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Abstract 


We  studied  the  capability  of  a  Fourier  spectrometer  system  to  remotely  detect  trace  gases 
in  localized  clouds,  e.g.,  stationary  source  effluents.  Detection  is  based  on  the  degree  to 
which  the  observed  IR  spectral  radiance  contrast  between  the  cloud  and  adjacent  background 
is  correlated  with  a  computed  reference  spectrum.  It  is  shown  that  trace  gases  can  be 
reliably  detected  when  spectral  features  are  well  below  the  noise  level.  The  minimum 
detectable  quantities  (MDQ’s)  for  various  trace  gases  at  one  atmosphere  total  pressure  are 
given.  The  MDQ’s  determine  the  combinations  of  gas  column  thickness  and  gas-background 
temperature  difference  that  correspond  to  95  percent  detection  probability  and  one  percent 
false  detection  probability  when  an  FTS  system  with  modeat-size  foreoptics  views  the  target 
through  a  path  equivalent  to  approximately  one  air  mass. 


Introduction  r 


Infrared-active  trace  gases  in  a  localized  cloud  can  be  remotely  detected  on  the  basis 
of  the  contrast  they  produce  in  an  IR  scene  image.  Of  course,  the  gas  temperature  and 
average  background  brightness  temperature  must  be  different,  and  the  resultant  spatial 
radiance  contrast  must  be  distinguishable  from  background  clutter.  Resolution  of  the  .image 
into  many  spectral  elements  will  give  additional  information  that  can  be  used  to  suppress 
the  effects  of  background  clutter,  atmospheric  spatial  variations  and  sy.'tem  noise,  and  to 
distinguish  between  different  trace  species  in  the  cloud. 

The  imaging  spectrometer  gives  the  spectrum  of  the  difference  between  target  (trace  gas) 
and  background  radiances,  which  can  be  compared  to  a  computed  reference  contrast  spectrum. 
Parameters  based  on  the  degree  of  correlation  between  the  two  spectra  can  be  uBed  to  decide 
whether  a  particular  target  species  is  present  and  also  to  characterize  its  amount. 

From  studies  reported  elsewhere1  <2  it  was  decided  that  a  mosaic  FTS  was  the  best  type  of 
sensor.  This  choice  was  based  on  an  assumed  target  cloud  size,  detection  range  and  sensi¬ 
tivity  requirements  (foreoptics  size  and  field  of  view),  the  number  and  spectral  locations 
of  the  target  species  IR  bands,  and  the  current  technology  of  high-performance  detector 
mosaics.  This  paper  gives,  for  a  particular  baseline  FTS  system,  the  minimum  detectable 
quantities  of  13  molecular  species.  System  noise-equivalent  spectral  radiances  for  the 
detection  bands  are  also  given,  which  allows  scaling  of  the  MDQ’s  with  respect  to  system 
radiometric  parameters.  The  reported  MDQ’s  correspond  to  idea)  conditions;  i.e.,  it  is 
assumed  that  system  performance  is  limited  by  detector  (system)  noise  rather  than 
atmospher xc/background  variations,  and  that  the  reference  spectrum  (which  includes  atmos¬ 
pheric  spectral  absorption  effects!  is  known  precisely.  The  effects  of  non-ideal  conditions 
are  discussed  in  Ref.  1. 


Theoretical  basis 

It  can  be  shown1  that  the  measured  spectral  radiance  contrast,  apart  from  system  spectral 
noise,  is 

ANV  =  N2  -Nj.  ~  DTvagv;  D  =  (Tg,Tsfc>  (1) 

where  tv  *  atmospheric  spectral  transmittance  between  the  target  and  sensor, 

agv  *=  spectral  absorption  coefficient  of  the  (single)  target  gas, 

u  «  molecular  column  thickness  of  the  target  gas  (number  of  molecules  per  unit 
area  in  the  line-of-sight) ,  and 
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Z¥v  -  difference  between  the  spectral  radiance  of  a  blackbody  at  the  gas  temperetu re 
Tg  and  the  actual  apectral  radiance  of  tha  background  at  tamperatura  T,fc , 
averaged  over  the  detection  apectral  bond. 

The  detectable  quantity,  denoted  by  D,  ia  essentially  a  acale  factor  in  the  measured  con¬ 
trast  spectrum.  Equation  (1)  is  an  approximation  only  because  we  have  replaced  the  actual 
spectral  variation  ABV  by  its  average  value.  However,  the  variation  is  ilow  compared  to 
that  of  rv  and  agv  and  is  nearly  linear  over  narrow  detection  bands,  so  that  the  approxima¬ 
tion  ia  a  good  one,  especially  if  tv  a  1  arid  agv  is  more-or-less  symmetrical  about  the 
center  of  the  detection  band. 

A  reasonably  accurate  theoretical  estimate  of  fvOgV  can  be  obtained  using  a  line-by-line 
computer  code  such  as  AFGL's  FASC0D1.3  Then  an  estimate  D'  of  D  can  be  obtained  by  finding 
the  value  that  minimizes  the  mean  square  difference  between  the  meaiurad  spectrum  ANV  and 
D*Tvagv  of  course,  ogv  is  a  function  of  the  gas  temperature,  but  31ra„v/3Tg  <<  SlnAB^/STg. 
Hence,  the  use  of  some  guess  temperature  in  computing  agv  will  nut  result  in  a  large  error’ 
in  D' . 


Some  simple  types  of  background  and  atmospheric  interferences  can  be  suppressed  by  sub¬ 
tracting  the  means  of  ANV  and  TvagV  from  these  spectra  before  they  are  used  to  estimate  D. 
The  least-squares  estimator  of  D  is  then 


£(a-a)  (b-E) 

£(b-&)2 


(2) 


where  a  =  ANV  +  N0ISEv ,  b  5  Tya_v,  and  the  summations  are  over  the  digitized  spectrum  values. 
Equation  (2)  is  the  prescription  for  computing  the  estimate  D'  given  the  measurement  a  and 
reference  ccntraat  spectrum  b.  By  combining  Eqs.  (1)  and  (2),  we  find  that  D'  can  also  be 
expressed  by 


u  «  +  — j 

% 

where  obn  is  the  covariance  of  the  digitized  reference  spectrum  b  and  the  spectrum  noise 
NOISE^ ,  and  ob2  is  the  variance  of  the  digitized  reference  spectrum. 

The  spectrum  noise  samples  will  be  Gaussian  with  zero  mean,  and  the  samples  will  also  be 
independent  provided  the  sample  spacing  is  approximately  equal  to  the  resolution  Av  of  the 
FTS  system.  It  can  be  shown  that  if  the  spacing  is  Av  the  second  term  of  Eq.  (3)  is 


cbn  x?n  /7  NESR 

-  « - *  x - =  xon,  (4) 

/Tz  rr:  _  U 

ob~  .M  ob  ob 

where  x  is  a  unit  normal  random  variate ,  M  is  the  total  number  of  resolved  spectral  ele¬ 
ments  in  the  spectra  a  and  b,  and  NESR  is  the  average  system  noise-equivalent  spectral 
radiance  over  the  detection  band.  Note  that  the  effective  rmB  noise  on  equals  /?  NESR 
because  ANV  is  the  difference  of  two  measurements  with  independent  noise.  Thus,  this  rms 
uncertainty  in  D'  is  oD-  5  /?  NESR/(/M  ob)  . 

The  quantity 


obD 

rT.  NESR 


.5) 


is  the  rms  signal  spectral  variation  over  the  rms  spectral  noise,  or  the  spectral  signal-to- 
noise.  Since  D/cD-  is  this  quantity  times  /H,  it  is  evident  that  for  sufficiently  large  M 
the  detectable  quantity  can  be  estimated  accurately  even  if  local  spectrum  features  are  wel1 
below  the  noise  level  (S/N  <<  1)  . 


D  will  be  negative  if  the  target  gas  is  observed  in  absorption  (is  "cooler"  than  the 
background)  and  positive  if  it  is  observed  in  emission.  Hence,  detection  of  the  target  gas 
can  be  based  on  the  test 


where  RD  is  some  predetermined  threshold.  Knowing  the  statistics  of  D'  (Eqs.  3  and  4)  we 
can  work  out  the  statistics  for  |D'|  for  both  target  present  (D  >*  0!  and  target  not  present 
(D  »  0)  ,  and  finally  determine  a  threshold  Rq  that  yields  satisfactory  detection  and  false 
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detection  probabilities .  It  can  be  ahown  that  tha  ttieofcold  value 

Rd  -  2.58  aD-  -  2.58  —  (7) 

D  D  /R  ob 

results  in  approximately  95  percent  detection  probability  and  approximately  one  percent 
false  detection  probability  when  D  equals • 

4.23  O  -  -  4.23  i  MOQ  .  (8) 

D  /H  ob 

This,  by  definition,  is  the  minimum  detectable  quantity. 

The  NESR,  M  and  ob  ail  vary  with  tie  spectrometer  resolution  Av.  A  near-optimum  value 
for  Av  is  the  halfwidth  of  the  target  gas  spectral  linen,  which  is  -0.1  crn'i  for  targets 
near  sea-level.  That  is,  Av  -  0.1  cm"l  will  give  the  smallest  MDQ  from  Eq.  (8).  The  single 
parameter  in  Eq.  (8)  that  depends  on  the  target  gas  and  detection  geometry  is  c^,  the 
standard  deviation  of  the  reference  spectrum  t degraded  to  0.1  cm'*  resolution. 

FTS  baseline  system 

The  FTS  baseline  configuration  that  was  selected  uses  i  four-port  Michelson  inter¬ 
ferometer  with  a  single  dichroic  and  two  detector  non, •Oca  in  each  of  the  two  output  ports 
(see  Fig.  1).  This  arrangement  allows  the  simulatneous  detection  of  up  to  four  different 
target  species.  The  fourth  port  (a  second  input  port)  provides  for  illumination  of  the 
backside  of  the  interferometer  beamsplitter  by  an  internal  blackbody  reference  or  by  a 
different  (background)  poition  of  the  external  scene  than  is  viewed  by  the  primary  input 
port.  Utilization  of  the  fourth  port  in  thin  manner  effects  an  optical  subtraction  of  the 
scene  and  reference  blackbody  spectra  or  of  the  displaced  scene  spectra,  and,  more  important, 
results  in  a  potentially  large  reduction  in  the  interfemogram  dynamic  range.'4  If  the 
reference  blackbody  and  background  have  the  same  brightness  temperature,  the  resultant 
target-background  spectral  contrast  obtained  by  this  optical  scheme  is  the  same  as  would  be 
obtained  using  a  conventional  single-port  interferometer  and  merely  subtracting  the  outputs 
of  two  different  mosaic  elements.  The  details  of  the  selected  baseline  configuration  are 
given  in  Ref.  1. 


MDQ  and  field  fit  results 
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The  MDQ's  were  calculated  for  several  candidate  detection  bands  of  each  of  14  different 
molecular  species.  The  "best"  band  for  each  species  was  selected  on  the  bav'xs  of  the  corres¬ 
ponding  minimum  detectable  column  thicknesses  (KIN  u’s)  for  a  given  small  temperature 
difference  between  target  cloud  and  background.  The  relationship  between  UN  u  and  MDQ  is, 
of  course,  strongly  dependent  on  the  spectral  location  of  the  detection  bead  through  the 
Planck  function  temperature  derivative.  Table  1  is  a  partial  list  of  the  beat  bands,  the 
NE§R’s  for  a  10-second  spectrum  scan  time,' the  corresponding  MDQ  and  the  MIN  v.  The  last 
quantity  is  the  minimum  detectable  concentration  in  ppmV  for  the  conditions  noted  in  the 
table  footnote. 


Gas 


HC1 

no2 

DF 

HBr 

n2o 

CO 

CH4 

S02 

NH3 

hno3 

C2H4°2 

ch3i 

NO-, 


_ Table  1.  MDQ  and  MIN  v  results _ 

,  *NESR  *MDQ  (molec/cm^) 

Band  (cm-1) _ (W/emr-sr-cm-^-) _ (W/cm^-sr-cm~i) 


2900-3000 

1.3 

X 

10-9 

2.5 

X 

109 

2650-2935 

1.4 

X 

10~9 

1.3 

X 

0 

t-i 

0 

H 

2700-2900 

2.6 

X 

10-9 

3.5 

X 

109 

2450-2650 

3.9 

X 

10~9 

2.4 

X 

1010 

2160-2210 

3.1 

X 

10~9 

3.1 

X 

0 

H 

O 

r”4 

2130-2185 

3.6 

X 

10*9 

7.2 

X 

1010 

1295-1310 

3.9 

X 

io~9 

2.1 

X 

1012 

1090-1210 

1.5 

X 

10“8 

2.7 

X 

1011 

915-970 

1.1 

X 

io's 

2.0 

X 

0 

H 

O 

H 

887.5-902.5 

6 . 0 

X 

10-9 

3.4 

X 

1010 

880-1000 

2.0 

X 

10*9 

5. 

X 

1010 

790-990 

2.4 

X 

10'8 

6.7 

:c 

1011 

710-795 

1.3 

X 

10'8 

4.0 

X 

1011 

**MIN  v  (ppmV) 
20 
92 
x7 
51 
22 
16 
165 
16 
1.0 
1.6 
2.8 
32 
19 


Note  that  the  MDQ  scales  as  the  system  NEST  (Eq.  8)  .  The  more  fundamental  quantity  ob 
for  each  detection  band  can  be  determined  from  Eq.  (8)  and  the  third  and  fourth  columns  of 
Table  1.  The  number  M  is  equal  to  the  width  of  the  band  divided  by  0.1  cm-*.  The  ob’s 
and  MDQ's  correspond  to  a  vertical  path  through  most  of  the  Hidlatitude  Summer  model 
atmosphere . 

Field  measurements  were  performed  at  Utah  State  University  to  val  ate  ths  theoretical 
predictions.  In  the  experiment  a  model  stack  emitting  a  controlled  now  of  NjO  in  front  of 
a  heated  "background"  plate  was  observed  by  an  FTS  system  at  a  distance  of  273  meters.  The 
detection  band  selected  for  the  experiment  is  2170-2260  cm' *  ,  which  includes  most  of  the  v-j 
band  of  N2O  (this  is  a  better  N2O  detection  band  than  the  one  given  in  Table  1  for  the 
measurement  geometry  of  the  experiment) .  The  NpO  waa  successfully  identified  and  quantified 
in  repeated  trials  when  th«  detectable  quantity  was  slightly  greater  than  the  calculated 
MDG.  Figure  2  shows  a  computed  (FASCOD1)  reference  contrast  spectrum  for  one  of  the  data 
runs,  and  two  examples  of  measured  contrast  spectra  that  resulted  in  N2O  detection.  The 
middle  panel  of  the  figure  shows  the  measured  contrast  when  the  detectable  amount  is  approxi¬ 
mately  four  times  the  MDQ;  in  the  bottom  panel  the  detectable  amount  is  -1.5  t  set  the  MDQ. 
There  is  extremely  high  probability  of  detection  (>>  0.95)  when  the  detectable  amount  is  1.5 
times  the  MDQ,  even  though  the  measured  contrast  spectrum  for  this  amount  looks  like  noise. 
The  spectrum  signal-to-noi se  defined  by  Eq.  (5)  is  0.26  for  the  bottom  panel  in  Fig.  2;  at 
the  MDQ  level  the  S/N  would  be  0.15. 
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Figure  2.  Normalized  reference  contrast  spectrum  for  V3  band  of  N20  (top  panel) 
and  two  normalized  contrast  spectra  obtained  from  measurements. 
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